Entanglement entropy and the Fermi surface 
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Free fermions with a finite Fermi surface are known to exhibit an anomalously large entanglement 
entropy. The leading contribution to the entanglement entropy of a region of linear size L in d 
spatial dimensions is S ~ L d_1 logL, a result that should be contrasted with the usual boundary 
law S ~ L 1 . This term depends only on the geometry of the Fermi surface and on the boundary 
of the region in question. I give an intuitive account of this anomalous scaling based on a low energy 
description of the Fermi surface as a collection of one dimensional gapless modes. Using this picture, 
I predict a violation of the boundary law in a number of other strongly correlated systems. 



Introduction. — Currently ideas from quantum informa- 
tion theory are stimulating major developments in quan- 
tum many body theory and vice versa. One of the ma- 
jor conceptual tools that has been receiving attention 
from both communities is entanglement entropy. To de- 
fine entanglement entropy the system must first be parti- 
tioned into two subsystems. The von Neumann entropy 
S = — Tr(pln/?) of the reduced density matrix of a sub- 
system is the entanglement entropy of the subsystem. 
If the quantum state of the system factorizes over the 
two subsystems then the entanglement entropy is zero. 
Also, if the whole system is in a pure state then the en- 
tanglement entropy is independent of which of the two 
subsystems is chosen. 

Bosonic systems frequently satisfy an area or bound- 
ary law: the entanglement entropy is proportional to the 
size of the boundary of the subsystem This bound- 
ary law can be interpreted as arising from short range 
entanglement. The boundary law is violated for confor- 
mal theories in one spatial dimension, and it is replaced 
by a logarithmic dependence S ~ logL on region size 
L (boundary law behavior would correspond to a con- 
stant independent of region size) @, H| ■ The coefficient 
of this term is universal and proportional to the central 
charge of the conformal field theory. However, conformal 
field theories in higher dimensions so far satisfy a strict 
boundary law. 

Interestingly, free fermions violate the boundary law 
even in more than one dimension A sharp d — 1 
dimensional Fermi surface coincides with an extra loga- 
rithmic scaling of the entanglement entropy [5|. Gapped 
fermionic systems and gapless Fermi systems with a 
higher codimension Fermi surface continue to obey a 
boundary law [1, 0] ■ A detailed formula for the entan- 
glement entropy of fermions with a sharp Fermi surface 
was proposed and found to be related to a conjecture of 
Widom important in signal processing [j| . In this Letter I 
give an intuitive picture of the anomalous entanglement 
entropy of free fermions with a Fermi surface in agree- 
ment with the Widom formula. The intuition is based 
on viewing free fermions as equivalent to a Fermi surface 
worth of 1 + 1 dimensional chiral modes. I also discuss 
other possibilities for observing violations of the bound- 



ary law. 

Free fermions and chiral modes. — I consider spinless 
fermions hopping on a square lattice in two dimensions. 
The lattice provides an ultraviolet regulator for the the- 
ory. It is possible to add spin and consider more generic 
lattices in different dimensions without any serious dif- 
ficulties. The physics is completely specified by giving 
the fermion dispersion relation e k as a function of the 
pseudomomentum k lying in the 1st Brillouin zone. In 
particular, the ground state is a Fermi sea where all states 
with energy less than the chemical potential fi are filled. 
For fermions at half filling there is one fermion per two 
lattice sites and hence the Fermi sea will occupy half the 
Brillouin zone. Generic filling fractions and dispersion re- 
lations lead to a Fermi surface with codimension 1. This 
means that the surface separating occupied and unoccu- 
pied regions in momentum space is a d — 1 dimensional 
subspace of the the d dimensional Brillouin zone. The 
generic presence of a finite density of states at the Fermi 
surface is responsible for many of the unusual properties 
of the Fermi gas. 

As an example of such an unusual property, the Fermi 
gas has a heat capacity linear in T in any dimension d. 
This result may be contrasted with the strongly dimen- 
sion dependent result for superfluid bosons C ~ T d . My 
interest here is in the anomalous entanglement entropy of 
the Fermi gas. Many bosonic and gapped fermionic sys- 
tems in d > 1 spatial dimensions have an entanglement 
entropy scaling as L d ~ x for a region of size L, but free 
fermions with a codimension 1 Fermi surface have the 
anomalous scaling L d_1 lnL. Remarkably, a precise for- 
mula for the logarithmic term in the entropy was conjec- 
tured and verified numerically 0^3- This formula reads 
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\n x ■ n k \dA x dA k , 
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where n x and n k are unit normals for the real space 
boundary and the Fermi surface respectively. The inte- 
grals are over a scaled version of the real space boundary 
(hence the overall L d_1 factor) and the Fermi surface, 
and the whole expression is written in units where the 
volume of the Fermi sea is one. I would like to under- 
stand how this anomaly arises from the presence of the 
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Fermi surface. 

Because of the finite density of states at the Fermi 
surface it is useful to reinterpret the low energy modes 
in terms of a large number of decoupled "radial" excita- 
tions. The word radial is used because the Fermi velocity 
Vk = Vfe£fc is normal to the Fermi surface. Thus I can as- 
sign to each point on the Fermi surface a radial fermionic 
mode with approximately linear dispersion and velocity 
given by the local Fermi velocity. These modes are ef- 
fectively relativistic and 1 + 1 dimensional traveling in 
the local radial direction as a function of time. They 
are also chiral because the direction of propagation is 
fixed by the local Fermi velocity. Modes traveling in the 
opposite direction are typically on the other side of the 
Fermi surface. Chirality is equivalent to the statement 
that there are no holes above the Fermi surface. This 
kind of patching procedure is the first step towards a 
higher dimensional analog of bosonization, but I will not 
need anything more than the heuristic picture of many 
chiral excitations. This point of view is also visible in 
various renormalization group treatments of Fermi liquid 
systems @-[l(|- 

Entanglement entropy. — The presence of a large num- 
ber of chiral one-dimensional excitations at low energy 
strongly suggests a violation of the boundary law. To 
make this precise I need some way to count the effec- 
tive number of such chiral modes. This counting can be 
performed using some intuition from the study of entan- 
glement entropy, namely that the imaginary partitions 
introduced to compute the entanglement entropy behave 
very much like real physical boundaries. Tracing out de- 
grees of freedom outside a region of characteristic size L 
should coarse grain the Fermi surface into patches of typ- 
ical size \/L d ~ 1 exactly as in a finite size system. Each 
patch contributes a factor of log L to the entanglement 
entropy because of the presence of a gapless one dimen- 
sional mode. The total contribution from the Fermi sur- 
face should therefore scale as L d ~ 1 \ogL with the Fermi 
momentum making up the units where necessary. This 
simple argument is made more precise below, but it al- 
ready captures the basic intuition behind the violation of 
the boundary law. 

Consider the case of half filling with dispersion efe = 
— 2eo cos k y a arising from a strongly anisotropic lattice 
model where fermions are unable to hop in the x direc- 
tion. The Fermi surface is nested and has two discon- 
nected components given by the lines k y = ir/2a and 
k y = — 7r/2a where a is the lattice spacing. Let the sub- 
system of interest be a box-like region of dimensions Lx L 
aligned with the x and y axes. The mode density on the 
Fermi surface is L/2ir and the length of the Fermi surface 
is 2ir/a + 2ir/a = An / a for a total of 2L/a modes. Each 
mode is chiral with left and right central charges given by 
cl = 1 and cr = 0. Left and right arc defined locally on 
each patch relative to Fermi velocity, but the important 
point is only that one of the central charges is one while 




FIG. 1: A sketch of the two model systems considered. Box 
A shows the Fermi sea (in gray) of the strongly anisotropic 
model. The relative orientation of the Fermi velocity and the 
chosen real space region (in black) is shown. Box B shows the 
Fermi sea (in gray) at half filling for a fermion hopping on 
a square lattice. Again, the relative orientation of the Fermi 
velocity and the chosen real space region (in black) is shown. 



the other is zero. The entanglement entropy of a one di- 
mensional conformal field theory on an interval of length 
L is Cl + Cr log (L/e) [3, E|. Putting everything together I 
find a total entropy of (L/a) log L/3. In order to compare 
with previous work I use units where the volume of the 
Fermi sea is one. This requirement is Air 2 / (2a 2 ) = 1 (the 
2 is for half filling) giving a = \[2tx. I find a total en- 
tropy of S = 3 ^2^ L log L in agreement with the Widom 
formula. 

Another simple situation is the case of equal hopping 
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with dispersion e k = — 2eocosfc x a — 2eocosfcj,a. At half 
filling the Fermi surface consists of a square rotated by 
45 degrees occupying half the Brillouin zone. To again 
make things simple consider an L x L region rotated by 45 
degrees so that it is oriented identically to the Fermi sur- 
face. Similar mode counting arguments now give 2\f2Lja 
modes. Using units where a = v27r, I find a total entropy 
of S = -r^L log L again in agreement with the Widom for- 
mula. 

After these examples an understanding of the general 
Widom formula can be obtained by breaking the real 
space boundary into small segments. I focus on the two 
dimensional case to make the notation as simple as possi- 
ble. Consider a segment AA X of the real space boundary 
and a segment AA k of the Fermi surface. With a mode 
density of AA x /2ir the patch AAk contributes 



AA x AAk\n x ■ n k \ 
2tt 



(2) 



modes. The flux factor \n x ■ n k \ counts the number of 
modes perpendicular to the real space boundary. In d 
dimensions the above formula is modified by replacing 
2ir with (27r) d_1 since AS X is now a general d— 1 dimen- 
sional surface element. Each of these modes is chiral and 
contributes 
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to the entanglement entropy on an interval of length L. 
Note that the precise choice of linear size L in the loga- 
rithm is not critical as differences can be absorbed into 
the non-universal boundary law piece of the entangle- 
ment entropy. 

The total entanglement entropy is given by integrat- 
ing the contributions from all patches of the real space 
boundary and the Fermi surface. This result must mul- 
tiplied by an additional factor of 1/2 because intervals 
are double counted in the integration. Equivalently, each 
end of a one dimensional interval can be thought of as 
contributing CL ^ R log (L/e) to the entropy 0j- The full 
result is 



S 



log L I 
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dA x dA k 
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and this is the Widom formula before rescaling the real 
space integral. 

Discussion. — I have shown in some simple cases that 
the entanglement entropy for free fermions can be ob- 
tained using counting arguments and intuition from one 
dimensional systems. I also argued that an explicit for- 
mula based on the Widom conjecture is the correct gen- 
eralization to arbitrary Fermi surface shape and region 
geometry. These arguments naturally reproduce the fully 
one dimensional case exactly. It is remarkable that the 
central charges of the chiral fermion are just right to make 
the Widom formula true. Nevertheless, the value of the 
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FIG. 2: Box A shows a spherical real space region of linear size 
L along with a chosen real space patch on the boundary. The 
size of the patch (bold black bar) has been exaggerated for 
clarity. Box B shows a local piece of the Fermi sea with filled 
states in gray. The dots are an exaggerated representation of 
the effective mode quantization coming from the real space 
patch. The relative orientation of the Fermi velocity oc n k 
and the real space patch normal n x is the origin of the "flux 
factor" \n x ■ n k \ counting the effective number of modes that 
propagate perpendicular to the real space patch normal n x . 
Here the angle between n x and n k is 7r/4. 



present point of view is not that it provides a rigorous 
derivation but that it gives intuition to search for other 
violations of the boundary law. 

The theory described here predicts a violation of the 
boundary law in any system with a large number of gap- 
less 1 + Id modes. Fermi liquids certainly qualify. The 
recently proposed d-wave Bose liquid phase should also 
violate the boundary law 



111 ]. Similarly, non- Fermi liq- 



uids with a sharp Fermi surface but no Landau quasipar- 
ticle should violate the boundary law [12j. Gauge/gravity 
duality has recently provided an interesting example of 
such a strongly correlated non- Fermi liquid phase [l3j . 
Frustrated quantum magnets where the low energy de- 
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scription is in terms of deconfined spinons with a Fermi 
surface should also violate the boundary law. In fact, 
this violation may be a useful numerical test for the ex- 
istence of such a spinon phase. This is because the log- 
arithmic correction to the boundary law is a low energy 
phenomenon that grows faster with region size than the 
non-universal boundary law term. 

Finally, I mention a few other implications of the the- 
ory. On the numerical side, the systems mentioned above 
should be difficult to simulate using tensor network ap- 
proaches 14 - TiJ . Such schemes are tailored to boundary 
law behavior in the entanglement entropy, and no nat- 
ural alternative is available for systems that violate the 
boundary law in more than one dimension. On the other 
hand, the violation of the boundary law could be a useful 
numerical signature of such phases using exact diagonal- 
ization or other methods that are not primarily limited 
by entanglement. 

Another implication concerns the appearance of the 
boundary law violating term across a phase transition 
from a gapped phase. The gapped phase near the phase 
transition can be viewed as a Fermi surface worth of 
gapped one dimensional modes. The simplest example 
would be a superconducting to Fermi liquid phase tran- 
sition where £ is related to the inverse of the supercon- 
ducting gap. Replacing L by the correlation length £, 
the gapped modes contribute an entropy U 1 ^ 1 log£ and 
obey the boundary law. As the quantum critical point is 
approached, the diverging boundary law term in the en- 
tropy is cutoff when £ exceeds L. The result is a bound- 
ary law violating term of the form L rf_1 logL as before. 
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